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The transition line of superconducting arrays of holes exhibits a rich field 
structure due to the interference of superconducting states nucleated at the 
holes edges. We studied by means of resistance measurements their effect 
on the T*(H) line as a function of transverse magnetic field using regular 
arrays of nano fabricated micron size holes. The arrays transition fields are 
higher than for the bulk. Moreover we found a nontrivial field modulation 
of the T* (H) line with an inversion, with increasing field, of the modulation 
concavity which we assigned to a crossover from a collective to an isolated 
edge state regime. The high field regime is well described by the nucleation at 
a single hole in an infinite film. The modulation at low fields was found to 
be dominated by the interference of neighbor edge states when the inter-hole 
distance w becomes comparable to the coherence length ^(T*). A comparison 
between arrays of different hole shape shows the influence of geometry on the 
type of interaction established, which can described either as a superconduct- 
ing wire network or as a weak link array. 

PA CS numbers: 74 . 76, 74 . 60, 74 . 25 



1. Introduction 



The effect of artificial pinning centers on the vortex dynamics of type- 
II superconductors has attracted a great interest in both fields of high T c 
and conventional low T c superconductors. Low T c materials are particularly 
interesting due to the development of nanofabrication techniques which en- 
able the introduction in a controlled way of nanofabricated defects. These 
patterned samples constitute relatively simple systems for the study of how 
geometric parameters such as defect size, shape and/or defect density influ- 
ence vortex dynamicsEJtrti 

Previous experiments on superconducting arrays of holes have shown 
that the transition line is dominated by surface superconducting states nu- 



cleated around the hole boundaries u. As a result, the array transition 
occurs at fields H* 3 higher than the nucleation fields in the bulk H C 2, and 
the transition line T*(H) exhibits a non-trivial field modulation due to flux 
quantization of the edge states over the two geometric lengths of the prob- 
lem: the hole surface and the array elementary cell. In particular, it was 
identified a crossover with decreasing field from an isolated hole regime to a 
low field regime dominated by the interaction between edge states. Though 
the behavior at high fields was well understood, the nature of the interactions 
at low fields remained unclear. 

In this paper we address in detail the collective regime at low fields 
when the distance between adjacent holes edges w is comparable to the 
array lattice constant a. The relevant parameters are the external field, the 
temperature dependent coherence length £(T*), the lattice constant a, the 
hole size, and the inter-hole distance w. We will try to put forward the role 
of these parameters on the crossover to the collective regime and the type 
of interaction established at low fields. The two extreme cases, the isolated 
hole in an infinite film! and the superconducting wire network §! (vanishing 
hole separation) are quite well known. 

In section ^] we present transport measurements on a superconducting 
square array where the aspect ratio hole size/lattice parameter is approx- 
imately two. These results are compared in section |3] with our previous 
work HI on a different array with similar aspect hole/lattice parameter but 
with twice the distance between hole edges. The distinct behavior found 
at low fields will be discussed as arising from distinct nucleation processes, 
determined by the ratio w/£. In section |3.2.| we analyze the case where 
w/^{T*) is lower than the critical ratio 1.84 0. Nucleation is then dominated 
by thin wire superconducting edge states and the coupling between them is 
well described by the formalism of superconducting wire networks II. In 
section |3.3. we analyze the case w/^(T*) > 1.84 where edge nucleation is de- 



termined at the single hole edge but w is still small enough to allow a weak 
overlap between neighbor edge state wave functions. In this regime the low 
field behavior can be described using a model based on a weak link interac- 
tion between adjacent edge states. This simple model is able to capture the 
essential features of the crossover previously discovered 0. 



2. Experimental Details 

The array sample consists of a thin film of aluminum (80 nm) patterned 
with a regular square array of nanofabricated holes. The lattice spacing is 
4.0 yum and holes have a square shape with 1.85 ± 0.01 fim side length, 
the distance between neighbor hole edges being w = 2.15 ± 0.01 /im as 
determined by SEM microscopy. We shall refer to this sample as sample A. 
The full array size is 1 x 1 cm 2 corresponding to a total of 6.25 x 10 6 holes. 

The patterning was defined on a monolayer PMMA coated 2" SLwafer 
by Deep UV photolithography using a high precision chromium maskB. The 
sample is then prepared using conventional lift-off techniques after thermal 
evaporation of pure aluminum over the resist mask on a UHV chamber. 
An homogeneous thin film of aluminum evaporated at the same time and 
submitted to the same fabrication steps is also measured for a reference of 
the patterning process effect on the material parameters. 

Both samples were studied by conventional four-probe resistance mea- 
surements using an AC four terminal resistance bridge at a 33 Hz frequency 
and a measuring current of 2 nA. Assuming an uniform current distribution 
over the array, the current density per wire is 4.5 x 10 Acm -2 . We used a 
similar in-line geometry for the voltage contacts (spiral-shaped) in the array 
sample as in our previous work H, placed at a distance of 2.8 mm from each 
other in the array center to avoid short circuits from the sample borders. For 
the reference sample we simply attached gold (non-superconducting) wires 
by ultrasonic bounding. No overshoot in the R(T) curves was observed, 
in contrast to the previously studied reference sample § where the spiral 
contacts were used and small cusps in the R(T) curves were present. [] 

Resistance measurements as a function of temperature were also per- 
formed at several magnetic fields between and 5 mT. From the zero field 
R(T) measurement we estimated several material parameters. Using a two- 

1 We remark that the R(T) curves for array A in non zero field do not exhibit the double 
transition observed on array B, which was then assigned to a first transition at T C 2 and 
to a second transition at higher temperatures T* . We now believe the observation of the 
T C 2 transition in array B was possible due to the transition of the 20 fim strips between 
patterned fields (of about 300 x 300 /im 2 ). In the case of array A, which is an homogeneous 
pattern over a 1 cm 2 surface, the double transition is not observed. 



dimensional Aslamazov-Larkin fit we estimate a BCS transition temper- 
ature T co = 1.262 K and a normal state resistance of R n = 0.099 f2. The 
resistance per single wire of length a = 4.0 /xm and width w = 2.15 /im is 
then r n = 0.354 and the normal state resistivity p n = 1.49 /ificm. Using 
vjp = 2.03 x 10 s cm/s we obtain an electronic mean free path l e i = 26.5 nm. 
The zero field transition temperature defined at half of the normal state 
resistance are 1.263 K and 1.265 K for the patterned sample and the ho- 
mogeneous thin film, respectively, the transition width being 3 mK for both 
samples. We summarize on Table 1 some parameters of array A and the 
parameters of the sample previously studied H, array B. This array consists 
on an aluminum thin film (80 nm thickness) patterned with a square array 
of circular holes of diameter 2r = 4.26 /xm and a lattice parameter of 9.0 p, 

The field dependence of the nucleation temperatures T*(H) is deter- 
mined using a heating feedback technique that keeps the sample resistance 
at a constant value while the magnetic induction field is smoothly varied by 
small increments of 0.2 /j,T. This method enables us to attain a fine field 
tuning of the array transition line. Several resistance criteria between 0.01 
and 0.7R n were used. The nucleation temperatures using the sweeping field 
method agree by less than 1 mK with those obtained for the same resistance 
criteria from R(T) measurements, indicating a good temperature regulation 
attained by the feedback method. 

The transition line T*(H) of array A for criteria 0AR n is displayed in 
Fig. |l|.a) as a function of applied field. For comparison it is also displayed the 
reconstructed bulk transition line, H c2 (T* /T co ) = (1 - T* /T co ) $ o /2vr£ 2 (0), 
using the reference film coherence length at T = K, £(0) = 220 nm (deter- 
mined from the initial linear slope of the reference film transition line). The 
array transition line is clearly above the estimation for the bulk. Besides, 
the non-trivial field modulation associated to the single hole and collective 
regimes are clearly identified. At fields below 0.75 mT it is characterized by 
periodic upward cusps, with a magnetic period of exactly H = 0.128 mT, 
that correspond to a quantum flux enclosed on a square cell of side length 
a = 4.0 //m, indicating the presence of phase coherence over the array lat- 
tice. At higher fields, the transition line exhibits downward dips occurring 
with a larger magnetic period. These large period oscillations correspond to 
the single hole regime discussed previously. In this case, the magnetic pe- 
riod is not constant since it depends on the effective surface formed by the 
hole radius and the surface superconducting sheath surrounding the hole. 
The average magnetic period is 0.543 mT, which corresponds to an effective 
square surface of side length Ihole = 1-94 /im, comparable to the array hole 
size of 1.85 fim. 

These field modulations were found for all the criteria used in the tran- 



sition line measurement. The crossover between the two field regimes is well 
illustrated on Fig. |l].b) where the sawtooth variation of the transition line 
slope dT* (H) I dH , at low field gives place to a smoother variation at higher 
fields, along with the change of magnetic period. 



3. Discussion 



3.1. Extraction of the energies for nucleation of 
superconductivity 

The role of the nucleation processes involved will be discussed in terms 
of the energies for nucleation of superconductivity at a given field, e nuc [. The 
energy e nuc i can be obtained by finding the lowest eigenvalue solution of the 
linear Ginzburg Landau differential equation 
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which fulfills the given boundary conditions on the order parameter ^(r). 

This approach is valid when we can neglect spatial variations of |vH r )|> 
such as thin films or wires of thickness <C £ or when the applied magnetic 
field reduces IVK 1 ")! to a value much smaller than the equilibrium amplitude 
\tpoo\ achieved deep inside the bulk superconductor. The regime of validity 
is then usually restricted to temperatures close to T co . 

An alternate approach is to extract e nuc i at a given field from the mea- 
sured T*(H) using the relation 00, 
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where ^> £>{x) = Y'{x)/Y{x) is the digamma function. This relation describes 
the depression of T*(H) relative to T co due to a magnetic perturbation. 
Though it was initially established by Abrikosov and Gor'kov for magnetic 
impurities 0, it can be generalized to all pair breakingnerturbations which 
destroy the time reversal symmetry of Cooper pairs EJ'EJ, if the scattering 
time of the electron pair is short enough for their relative phases be random- 
ized by the perturbation. 

It can thus be applied to a dirty superconductor in strong external 
magnetic fields (and only surrounded by insulators), if the mean free path 
l e l is much smaller than all sample dimensions and £ D or in the case of a small 
superconducting particle with all dimensions <ti £ a . The sample thickness 
must be smaller than £ and the penetration depth to ensure the penetration 



of the magnetic field. At temperatures close to T*(H), e nuc i is the energy 
required to break the Cooper pair thus destroying superconductivity. When 
T* = 0, (or H = iJ* 3 (0)), € n uci coincides with the BCS superconducting gap 
1.76kBT co . Close to T co the digamma function can be expanded around 1/2 
and the T*(H) depression is linear in e nuc i, T co — T*(H) = e nuc iir /8ks- 
The advantage of using Eq. ^ on the determination of a n uci is 

that it 

remains valid down to all temperatures and in strong magnetic fields. The 
linear Ginzburg-Landau results can be recovered if the temperature depen- 
dent coherence length is defined as £ 2 (T) = Dfo/e nuc i, where D = l/3vpl e i 
is the coefficient for electronic diffusion. For sample A, D = 180 cm 2 s ob- 
tained from the mean free path l e [ = 26.5 nm. Using £ 2 (0) = hD/1.7QkBT co 
we estimate £(0) = 250 nm. 

We calculated e nuc i(H) for arrays A and B using Eq. || and the experi- 
mental T* at the given applied field H. These results are displayed in Fig. § 
after being normalized by the nucleation energy on the bulk, e C 2 = hDH/<& , 
with the same coherence length as the array. For comparison, we also rep- 
resent the theoretical e nuc ;/e C 2 for a circular hole on an infinite film (solid 
line) 111. In fact, the representation e nuc i/e C 2 is equivalent to the inverse ra- 
tio of the nucleation fields H C 2/H* 3 , that close to T co acquires the Ginzburg 
Landau form H c2 /H* 3 = {l-T*/T co ) H c2 {0)/H, using H c2 {0) = $ o /2vr£ 2 (0) 
and the coherence length as defined above. 

From Fig. |2] it is clear that in the high field regime both samples are 
very well described by the theoretical single hole case. With decreasing fields 
both arrays deviate from the single hole description, with the appearance 
of the collective field modulation, periodic on <3? per array cell, and with 
upward concavity. It is in this regime that the samples present a strikingly 
distinct behavior. The reduced energies e nuc i/e C 2 for sample B and for the 
single hole tend to an overall increase with decreasing field, reaching 1 at 
zero field. In contrast, for sample A, e nuc i/e C 2 decrease with decreasing field, 
dropping well below the single hole line. This means that in this regime, 
array A presents a ratio H*a/H C 2 higher than the classical limit of 1.69 
for an infinite surface sheathQ. Comparing the array geometric parameters, 
they both have similar aspect ratios w/a, 0.54 (array A) and 0.53 (array B), 
respectively and similar ratios of the superconducting volume over the array 
cell volume V s /V ce u, 0.79 (array A) and 0.82 (array B). We thus believe the 
distinct low field behavior is associated to the different ratio w/^(T*), which 
close to T co controls the process of edge nucleation and the type of coupling. 
In fact, both samples loose the collective behavior for w/£(T*) > 3. 

On the following subsections we shall analyze these features by describ- 
ing the array nucleation energy as coming from two main contributions: the 
nucleation energy of the single edge state and the coupling energy between 



neighbor edge states. 



3.2. Wire networks of wide strands 

In this section we focus on the low field behavior of array A. In this 
regime, the array energy is a sum of the nucleation energy on a individual 
wire in parallel field e strip \\ and a coupling energy, that can be described 
within the framework of superconducting wire networks theory BE. 

The case of superconducting wire networks of narrow wires is well un- 
derstood. It can be treated as a periodic array of superconducting islands 
strongly coupled to each other by thin superconducting wires of length a 
(£ » a) and width w « £. Neglecting superconducting fluctuations, the 
coupling energy can be computed within mean-field theory by solving the 
linearized Ginzburg-Landau equations at each node of the network. For a 
regular square lattice (same length of all strands) the order parameter tpi 
on each island i will be coupled to the first neighbor sites through field 
dependent phase factors as, 

4 cos (u) ipi = ^2 4>j exp(-i7 ii ) (3) 
5 

where 4 is the lattice coordination, 7^ = 2tt/& f? A • d\ is the phase factor 
along the wire linking site i to a site j, A the vector potential, and u = 
a/£(T) is the strand length in units of The coupling energy is then 

given by, e wnt = TiDu 2 /a 2 . 

Eq. H] is equivalent to a tight binding equation in a potential with the 
same geometry and tight binding energy etb = 4 cos u. The network coupling 
energy e wn t = HDu 2 /a 2 can then be expressed in terms of the tight binding 
ground state level etb with u = arccos (e^/4) 00. However, this relation 
is only valid in the limit w >> £. Taking into account the finite thickness 
of the wires, a more complex relation between etb and u is obtained, etb = 
4 cos u + 4 t&n(uw/2a) sinu — uw/a 113. This result was established for 
zero external field but it can still be applied at low fields while there are no 
vortices in the wires. 

In the inset of Fig. || is displayed the coupling energy of array A and the 
theoretical e wn t energy as a function of magnetic flux per array elementary 
cell <3? ce /z, for & ce u/<& between and 1. The coupling energy for array A 
is obtained from e nuc i (calculated using Eq. [|| and the experimental T*(H)) 
after subtracting the parabolic energy contribution due to edge nucleation 
in the wires, e strip 11 . The theoretical e wn t is obtained from the ground state 
tight binding energy calculated for rational values of Q C eii/&o = p/q, with 



q < 30 and p < q. A very good agreement is obtained between the experi- 
mental and the theoretical data which take into account the wire thickness. 
The theoretical results in the limit w = lead to a smaller energy. Since 
we are close to T co , the variation of T co — T*(H) due to the coupling (after 
subtracting the linear dependence on H) corresponds to e wn t7r/8kB- 

Besides the fundamental dips at & C eu/^o = 0, 1 and at 1/2, additional 
dips can be identified at rational values $> ce u/& = p/q, for q = 3,4 and 
5. This field structure is a manifestation of the interference of quantum 
states over cells of size qa x qa 111. With increasing field, the fine field 
structure becomes less pronounced and only the fundamental dips remain 
until &ceii/&o = 8. At higher fields the single hole regime is recovered. 

In fact, the crossover to the single hole regime is associated to a crossover 
from a two-boundary to a one-boundary nucleation process in the network 
wires. This explains the smaller values of e nuc i/e C 2 when compared tothe 
classical limit e nuc z/ e c2 = 0.59 for nucleation in an infinite surface sheatrJH or 
for nucleation in a single hole. Such as for a thin slab in a parallel field, edge 
nucleation in the array at low fields is controlled by two boundary conditions 
imposed at the edges of adjacent holes. 

The field dependence of e str i p y for a strip or slab of intermediate thick- 
ness w, is strongly dependent on w/£(T*). Below a critical thickness w < 
1.84 £(T C *) (thin wire limit) nucleation starts symmetrically at both surfaces 
and the maximum of the order parameter occurs at middle distance between 
them. In this limit e str i p y = H 2 w 2 (irhD/6^). With increasing field, when 
w > 1.84 £(T*), e strip || deviates from the parabolic field dependence as 
the order parameter solutions at each surface pull apart, their sobreposition 
giving rise to nodes along the middle plane of the wire and equidistant of 
AL ~ <& Q /Hw(l — 1.84£). With increasing field the vortex pattern becomes 
more complex, until the interference between the surface solutions become 
negligible (compared to ksT*) and the one-boundary solution is recovered. 



t 0.59hDH/$ o . 
Here we are inter- 



In this limit e str i p y approaches the surface sheath resu : 
All these features were discussed previously 00111 . 
ested in comparing the envelope of the energy curve for array A and the 
nucleation energy for a strip, e strip m. On Fig. || is displayed e nuc i for ar- 
ray A and e stT i p ii for a wire with the same width w as the array strands 
as a function of applied field. We can identify several similitudes. On the 
single hole regime the main dependence of e nuc i is linear on H such as the 
field dependence of e str i p y in the surface sheath limit, in agreement with the 
dominance of one-boundary nucleation at the edges of each individual hole. 
With decreasing field both curves deviate from the linear field dependence 
due to the emergence of interference between adjacent surfaces solutions. 
This deviation occurs near the field H\ = 0.65 mT which corresponds to the 



position of the first important collective dip of e nuc i/e C 2 at Q ce u/& = 5 (see 
also Fig. £§b). Below the field H a 2.75 & /irw 2 (0.39 mT) nucleation starts 
symmetrically^ in the strip and the order parameter is maximum at middle 
distance between adjacent edges. The occurrence of symmetric nucleation 
on a low field regime explains why the wire network description is still valid 
for arrays of wide strands. In the regime H < H Q , the collective dips of e nuc i 
approach the parabolic envelope e strip m, since at integers values of <& ce u/<& 
the costs in coupling energy are minimum. These results thus indicate that 
the small values of e n ud/ e c2 for array A at low fields are simply related with 
the two-boundary nucleation process. 

On the other hand, since the crossover from two-boundary to one- 
boundary nucleation is associated to the appearance of interstitial nodes 
of the order parameter within strands, we expect a broadening of the array 
resistive transition with increasing field due to these weakly bounded in- 
terstitial vortices. In Fig. |4] is represented the field variation of the resistive 
transition width AT*(H), for array A obtained by subtracting the transition 
lines T*(H) measured for resistive criteria 0.7R n and 0.03i? n . For compari- 
son, we also represent the field variation of the distance AL between nodes 
on a single wire, normalized by the array lattice parameter a = 4 (im. 

At fields H < H = 0.39 mT (region I, <5> ceU /<S> < 3) the transition 
width at integers & ce u/& is the same as in zero field since there are no 
order parameter nodes in the strands, only coreless vortices inside holes. At 
integers <& ce u/& , the flux quanta per array cell corresponds to the quanta 
enclosed by each hole. The transition is then sharpened since every hole 
encloses the same number of flux quanta. At intermediate values of & ce u/& 
phase fluctuations lead to a small broadening of the transition (of about 
1 mK). 

At fields H a < H < Hi (region II, 3 < § ce ii/<&o < 5), the first nodes of 
the order parameter are expected to appear within the strands with a sepa- 
ration AL that drops from infinity to values comparable to the array lattice 
parameter. In this regime there will be a competition between increasing the 
flux per hole by <I> or follow the increase of field by accommodating vortices 
at interstitial positions in the array wires. The increase of the transition 
width for $ ce ///<I> between 4 and 5 illustrates the presence of a few loosely 
bound vortices. This observation is in agreement with Fig. |].b) where the 
field dependence of e nuc i/e C 2 represented as a function of & ce u/Q shows that 
e n ucz/e C 2 at ® ce u/$o = 4 is higher than the adjacent dips. At § ce u/§o = 5 
every hole encloses 4 <£ which favors the decrease of e nuc //e C 2. Also, the 
distance between vortices in the strands should be of the order of a and, in 
that case, they can occupy positions at the array interstices forming a stable 
sub-lattice, which decreases the resistive transition width. 



In region III, the distance between interstitial vortices should drop be- 
low a and the transition width broadens considerably due to these weakly 
bounded vortices. The distance between nodes is further reduced with in- 
creasing field until edge states at each hole become independent. 

We can see from Fig. |].b) that at & C eii = 8 $ 0) the field modulation of 
e nucl/ e c2 inverts concavity once the single hole nucleation becomes dominant. 
If for $ce^/$ between 5 and 8, the additional flux occupies interstitial 
positions, the flux per hole remaining equal to 4 $ D , e nuc [/e C 2 should meet 
the single hole curve at its A th oscillation, in agreement with the results 
presented on Fig. |2[ 



3.3. Weak link array of edge states 

We do not expect the description of a wire network of wide strands to 
hold if holes edges are pulled apart and/or shaped to circles. This is the case 
of sample B where the minimum inter-hole distance is w = 4.74 fim ( ~ twice 
compared to array A). The field for symmetric nucleation is H Q ~ 0.08 mT 
for a stripe with the same thickness. However since the distance between hole 
edges is not constant (varying from 4.74 to 9 //m around the hole perimeter) 
H should be further reduced. As a consequence the parabolic envelope of 
t nuc i is not observed. The coupling between neighbor edge states is still 
present though. Fundamental dips of e nuc i at integers values of <& C eii/^o < 7 
are clearly observed although we find no fine field structure £3. At low field 
e nuc i/e C 2 is lower than the single hole calculation indicating that the energy 
contribution due to the overlap between edge states is important. 

In the next sections we present a simple model which is able to capture 
the transition from the single hole to the collective behavior for array B. The 
application of this model is based on the assumption that in a narrow tem- 
perature region T C 2 < T < T* the superconducting edge states are localized 
close to the hole boundaries and the system behaves like an array of weak 
links. The fast oscillations observed at low fields are then related to the 
collective behavior of the array. Since there is nonzero overlap between the 
wave functions of different edges, the energy should be higher than that of 
a single hole. One consequence is that the nucleation field H* 3 is lower than 
that of a single hole, as it is observed in the experiments. The single hole 
behavior is naturally recovered at higher fields where the overlap between 
neighboring edges vanishes. 

A Ginzburg-Landau approach is suitable to study the problem. For the 
single hole, an analytical solution is available due to the cylindric symmetry 
of the problem c3. A variational approach for the determination of the order 



parameter has been analyzed by Buzdin □. The agreement with the exact 
solution is good except at low fields where the variational approach predicts 
a region with zero flux through the hole which is not present in the exact 
solution. Buzdin's variational approach, however, has the great advantage 
of being applicable also in the many holes case where an analytical solution 
is impossible to obtain. 

We proceed as follows. We first describe an improved variational ansatz 
for the single hole case by using a trial order parameter of three parameters. 
We then construct an effective Ginzburg-Landau free energy which takes 
into account, in an approximate way, the edge coupling and we determine 
the resulting nucleation energy. 



3.4. Variational approach : single hole 

In this section we introduce the new variational wave function for the 
evaluation of the nucleation energy in a superconducting film containing a 
hole. Following Buzdin □, the Ginzburg-Landau free energy ( Fn is the free 
energy of the normal state) is given by 
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The fourth order term has be ignored since we are interested in the phase 
boundary. The surface critical field can be determined with a variational 
procedure by determining the minimum of the functional 

F-F N = (5) 

within the class of trial wave functions ip(f) which satisfy the proper bound- 
ary condition at the superconductor-insulator interface 

-dp\ P =R 

where R is the radius of the hole. For a hole containing m flux quantum, the 
order parameter has the following form in cylindric coordinates (p, 4>, z), 
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. It is possible to improve the result obtained by Ell by using a three- 
parameters trial order parameter of the form 



F(x) = 1 + a(x — xo) v + P(x — xq 
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where a, (3, 7 are the variational parameters and the dimensionless quantities 
x = p(eHfch) 1 ^ 2 , x = R(eH '/ 'cK) 1 / 2 = (^hoie/^o) 1 ^ 2 has been introduced. 

For the case of a single hole an exact solution for the nucleation energy 
is expressed in terms of the eigenvalues of the Kummer's equation We 
compare our results (the best results have been obtained for r] = 2 and £ = 
2.05) with the exact calculation and with the one parameter trial function □. 
This is shown in Fig. |5[ The trial function F(x) given by Eq. (|7|) has a 
maximum at (x — xq) 2 ~ (2 — 7/(0 + /?))/7, a feature which is also present 
in the exact solution. This seems to improve considerably the accuracy of 
the variational approach E3. 

The result obtained by using Eq. (|7|) improves considerably the agree- 
ment with the exact solution. In particular the spurious m = is shrunk to 
very low fields (< 0.03 $ ). 



3.5. Variational approach : array of holes 

We obtain an approximate expression of the nucleation energy for a 
regular array using the variational wave function introduced in the previous 
section and an appropriate Ginzburg-Landau free energy which takes into 
account the overlap between edge states of adjacent holes. The model is 
sketched in Fig. |(| we will consider an array of edges connected by weak 
links (dashed lines in the figure). The crossover between the single hole and 
the collective behavior stems from the interplay of flux quantization on each 
hole and frustration effects through the elementary array cell. 

The GL free energy which we propose has the following form 

•^-•^v = E^ + E^ ( 8 ) 

where T% is the free energy defined in Eq. (||) and the subscript i identifies 
the hole in the array. Tij takes into account the overlap between the holes. 
We choose a J-™ of the following form 
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a Q \ifii(f) - ^-(f)| 2 + (ai + a 2 5 ij ) |^(r) || Vj(r)|] (9) 



The effect of the external magnetic field will be included in Eq. |9] by means 
of a Peierls substitution. We assume that the overlap is restricted to a small 
area (see the dashed lines in Fig. ^) and that the hole array can be treated 
as a weak link array. The first term in Eq. ^ is modified as follows 

\ifii - ipj\ 2 -> IV'il 2 + \ipi\ 2 - 2\<ipi\\ipj\cos((f>i - (f)j -irm- Aij) (10) 



where A^j = (27r)/($ G ) f? A • di is calculated along the path indicated by 
the dashed lines in Fig. |6| (we assumed that all the holes contain the same 
number m of flux quantum). 

The terms in Eq. || can be understood considering that the order pa- 
rameter in the array of holes can be given as tp = J2i V'i • Substituting this 
expression in the GL free energy of Eq. |], one generates the various terms 
given above. In general since the approach is phenomenological, the various 
contributions enter with different coefficients (a Q , 01,02)- It is important to 
stress that the total wave function does not satisfy the proper boundary con- 
ditions around the holes because of the exponential tails of the edge states of 
the neighboring island. Therefore this approach breaks down at very small 
field when the overlap becomes too strong. 

The agreement of the experiments with the single hole result in the 
high flux regime indicates that phase fluctuations do not drive the phase 
transition. In this case we can approximate the phase dependent part of the 
free energy by its ground state energy 

2 

^ceii/^o) = ■^^2( cos ( ( t>i ~ <Pj - vrm - Aij)) GS 

where N is the number of holes in the array, z is the coordination number, 
^ceii/^o is the magnetic flux per elementary cell of the array and (. . .)gs 
means the ground state configuration of the .phases fa. 

Going over the same steps as Buzdin □ the nucleation energy for an 
array of holes is obtained by minimizing the following functional over the 
trial function introduced in Eq. ^ 



^ = zf° — 



Cc2 Jxq 



(m - x 2 fF'\x) + x 2 [F\x)} 2 + g ^-L x *[F{x)f 



+ 9iK^^-92e^> cell /^ )\{F{x i )F{x j )) (11) 

CI l -±Jbr\ 



where the overlap integral 



/w 
dxF(x)F(x + —) (12) 
x 



the normalization 



POO 

1 L = 2 / dx xF 2 (x) 

Jxo 



has been introduced. The new parameters <?,<?i,<?2 are easily expressed as a 
function of the a's. 



In Fig. ^ we present the results obtained for various values of the three 
parameters g,gi and g<i- We can observe that the presence of the Joseph- 
son coupling increases the array nucleation energy, relatively to the isolated 
hole case, and introduces cusps of different concavity superimposed over the 
main single hole background at low flux when the coupling becomes stronger. 
These features were found in the experiment of Bezryadin et al a. At very 
low fields e nuc i/e C 2 rapidly increases, as also seen in the experiments. It may 
be explained by the fact the the wave functions extends over various lattice 
constants and therefore the system does not show any surface superconduc- 
tivity. However this model breaks down at flux typically of the order of 
0.15 $ ho i e . 

A word of caution is needed at this point. Although the main features 
of the crossover are found one should be aware that the model is still too 
simplified to aim at a quantitative comparison with the experiments. In 
particular the choice of the values of the constants g is not related to any 
microscopic model which allow to justify the numerical parameters. It is, 
however, rewarding that most of the qualitative features are captured by 
our model. Afurther step might be to apply the techniques developed by 
Palacios et al Ej to the present problem. 



4. Conclusion 

In conclusion, we tried to show how geometric parameters such as the 
inter-hole and the array lattice constant can influence the transition temper- 
ature of periodic superconducting arrays in a magnetic field. 

We found a clear crossover with decreasing field from a single to a 
coupled edge state regime for both samples. The behavior of e nuc i/e C 2 in the 
single edge regime is very similar and the arrays transition depend exclusively 
on the magnetic flux per hole area. However, we remark that a quantitative 
comparison of the energy ratios e nuc i/e C 2 between samples need some caution 
since they are extremely dependent on the estimation of £(0) and small errors 
can lead to a shift of the e nuc i/e C 2 values. 

The low field behavior of the studied samples is representative of two 
distinct coupled systems: the superconducting wire network (strong cou- 
pling) and the weak link array. In both cases the coupling is mediated by 
the frustration induced by the applied field over the array elementary cell, 
however the type of coupling depends on the inter-hole distance w, the defect 
shape and the coherence length at T*(H) which controls the edge nucleation 
process. In the case of the weak link array (array B) we have a coupling of 
single hole edge states which depend on H through the enclosed flux per hole. 



However, in the low field regime, the enclosed flux does not necessarily follow 
the field increase since the change of the hole winding number may be less 
favorable than placing flux at the array interstices where superconductivity 
is weakened. In the wire network case (array A at low fields), we have a 
strong coupling of the order parameter at the array nodes, which is non-zero 
along the wires due to symmetric nucleation. The flux enclosed per hole or 
per array cell are equivalent quantities as long as interstitial vortices are not 
allowed within the array strands. In this regime, the main contribution to 
the array transition comes from nucleation in the wires (which depends only 
on H). 

This explains why the studied samples who have similar geometric ratios 
(w/a 2 and V s /V ce u « 0.8) exhibit such a different collective behavior. 
The relevant parameters are then the inter-hole distance (which allows two- 
boundary nucleation in a wide field regime for sample A) and the parallel 
hole edges which favor the thin wire nucleation. For periodic arrays of close 
enough (w/^(T*) < 1.84) and parallel hole edges, the coupling between edge 
states can be described as a wire network coupling and the array transition 
approaches the thin wire nucleation at low fields. As a consequence, the ratio 
H c z/H C 2 (e nuc i/e C 2) in these arrays is not necessarily upper (lower) bounded 
by the surface sheath result H C ^/H C 2 = 1.69 (e nuc //e C 2 = 0.59) or the single 
hole nucleation limit. 

When two-boundary nucleation does not occur, due to a higher inter- 
hole distance and/or due to the hole shape, the wire network formalism is not 
suitable to describe the coupling between edge states, since the variations 
of the order parameter amplitude between holes cannot be neglected. In 
this case the array transition is dominated by single hole nucleation and we 
presented a coupling description, based in a weak link interaction between 
single edge states, which is able to recover the main features of the T*(H) 
line on the low field regime: the inversion of field modulation concavity 
and the presence of periodic upward cusps. In addition, it reproduces the 
experimental observation that the array transition occurs at temperatures 
below the single hole nucleation, since the array energy is increased due to 
the overlap between neighbor edge wave functions. 
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Sample 


hole shape 


a (//m) 


w (/mi) 


T co [K] 


m (Mm) 


array A 


square 


4.0 


2.15 


1.263 


0.25 


array B 


circle 


9.0 


4.74 


1.25 


0.25 



Table 1. Some parameters of array A and array B. The array lattice constant 
is a and w the minimum distance between the edges of adjacent holes. 



Fig. 1. a) Field dependence of the superconducting transition line T*(H) of 
array A (solid line) and of the reference sample T C 2(H) (dashed line). Two 
types of field modulation are clearly identified for array A: downward, large 
period oscillations with dips at half integers of <£ per hole (down arrows) and 
upward oscillations of shorter period with cusps at integers of <J> per array 
elementary cell (inset :up arrows); b) T*(H) slope for array A as a function 
of H. The change of magnetic period due to the crossover from collective to 
single hole regime is quite visible. 



Fig. 2. Normalized nucleation energies e n ud/ e c2 a s a function of magnetic 
field (in units HShoie/^o = &hoie/&o), f° r sample A (solid dots), sample B 
(open diamonds) and the theoretical calculation for a cylindric cavity in an 
infinite thin film (solid line). 




H (mT) 

Fig. 3. Field dependence of the nucleation energies of array A (solid dots) 
and for a strip of width w = 2.15/xm (open diamonds), normalized by kBT co . 
In the field range H a < H < H\ interstitial vortices appear within strands. 
Inset: Coupling energy e wn t for array A (solid dots) and the theoretical e wn t 
for a superconducting wire network with w = (small dots; lowest curve) 
and taking into account the wire thickness (small dots, upper curve) as a 
function of reduced field <& ce u/& between and 1. The main dips position 
at rationals pjq are indicated by down arrows. 




Fig. 4. a) Resistive transition width AT* of array A (open diamonds) as 
a function of the reduced flux $ ce ^/$ , and comparison with the normal- 
ized distance AL/a between interstitial vortices for a thin wire (solid line). 
An oversimplified picture of the vortex patterns developed within the wires 
is represented. Three main regions can be identified: (I) w < 1.84£(T), 
nucleation starts symmetrically and there is no vortex in the wires; (II) 
w > 1.84£(T) and oo > AL/a < 1, nodes of the order parameter appear 
at interstices due to the interference of neighbor edge wave functions (white 
dots); (III) AL/a < 1 and decreases with increasing field until the surface 
solutions become independent and the single edge states are localized around 
each hole, b) field variation of the nucleation energy e nuc i/e C 2 for array A. 
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0.4 0.8 1.2 1.6 

^hole' ^0 

Fig. 5. Variational approximations for e nuc i/e C 2- The thick curve is ob- 
tained using the three parameter variational wave function. For comparison 
the curve obtained from the one parameter variational function is reported 
(dashed line). 




Fig. 6. Elementary cell of a square array of holes (grey circles). The edge 
states localized around each hole (thick dashed lines) are weakly coupled to 
first neighbor edge states. The coupling is indicated with crosses. 
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Fig. 7. The normalized nucleation energy e nuc i/e C 2 for a regular array of 
holes obtained by minimizing the functional as defined in the text for various 
parameter values: a) g=0.1 gl=0.4 g2=1.6; b) g=0.2 gl=0.7 g2=0.9. 



